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Abstract 

We present the results of an investigation to determine the longitudinal (zonal) distributions 
and temporal evolution of ices on the surface of Triton. Between 2002 and 2014, we obtained 63 
nights of near-infrared (0.67-2.55 pm) spectra using the SpeX instrument at NASA’s Infrared 
Telescope Facility (IRTF). Triton has spectral features in this wavelength region from N 2 , 
CO, CH4, CO 2 , and H 2 O. Absorption features of ethane (C 2 H 6 ) and are coincident at 

2.405 pm, a feature that we detect in our spectra. We calculated the integrated band area 
(or fractional band depth in the case of H 2 O) in each nightly average spectrum, constructed 
longitudinal distributions, and quantihed temporal evolution for each of the chosen absorption 
bands. The volatile ices (N 2 , CO, CH4) show signihcant variability over one Triton rotation and 
have well-constrained longitudes of peak absorption. The non-volatile ices (CO 2 , H 2 O) show 
poorly-constrained peak longitudes and little variability. The longitudinal distribution of the 
2.405 pm band shows little variability over one Triton rotation and is 97±44° and 92±44° out 
of phase with the 1.58 pm and 2.35 pm CO bands, respectively. This evidence indicates that 
the 2.405 pm band is due to absorption from non-volatile ethane. CH4 absorption increased 
over the period of the observations while absorption from all other ices showed no statistically 
signifcant change. We conclude from these results that the southern latitudes of Triton are 
currently dominated by non-volatile ices and as the sub-solar latitude migrates northwards, a 
larger quantity of volatile ice is coming into view. 
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1. Introduction 

The images of Triton taken during the Voyager 2 flyby in August 1989 provided a tantalizing 
view of a diverse and dynamic surface. Two adjacent terrrains on the imaged hemisphere 
appear drastically different in both color and texture, hinting that these two distinct surface 
units may have different surface ice compositions (Stone and Miner, 1989). Triton is frequently 


1 



compared to Pluto, and with good reason. Triton and Pluto are comparable in size and may 
have come from the same initial population; it is believed that Triton is a captured Kuiper Belt 
Object (Agnor and Hamilton, 2006). However, Triton has a more diverse collection of surface 
ices, with N 2 , CO, CH 4 , CO 2 , and H 2 O ices dehnitively identihed in ground-based spectra 
(Cruikshank et ah, 1993; Cruikshank et ah, 2000). Spectral signatures, possibly due to C 2 H 6 
(ethane), are present as well (DeMeo et ah, 2010). Two independent investigations suggest that 
Triton’s surface temperature is about 38 K (Broadfoot et ah, 1989; Tryka et al, 1993). At this 
temperature, the sublimation pressures of N 2 (22 /ibar), CO (3 /ibar), and CH 4 (0.002 pbar) 
are non-negligible, while the sublimation pressures of CO 2 ( 10 “^^ /ibar), H 2 O (0 /ibar), and 
ethane (2x10“^® /ibar) are negligible (Fray and Schmitt, 2009). Henceforth, we will refer to 
N 2 , CO, and CH 4 as volatile ices due to their relatively higher sublimation pressures compared 
to the non-volatile ices: CO 2 , H 2 O, and ethane. CO 2 and H 2 O are less mobile and presumably 
constitute the substrate upon which the volatile ices (N 2 , CO, CH 4 ) deposit. The presence of 
volatile ices on the surface, though predicted for an object with the surface temperature and 
diameter of Triton (Schaller and Brown, 2007; Johnson et ah, 2015), is surprising. The capture 
process and subsequent circularization of Triton’s orbit through tidal interactions heated Triton 
signihcantly, resulting in the production of a thick atmosphere and blowoff of volatile species 
(McKinnon et ah, 1995). Present-day surface composition is puzzling in light of the dynamical 
history of Triton. 

The continued presence of volatile ices allows for an atmosphere around Triton. Despite its 
low surface temperature and high geometric albedo (0.719; Hicks and Buratti, 2004), the surface 
pressure on Triton was measured at 14±1 /rbar by Voyager 2 in 1989 (Gurrola, 1995). Stellar 
occultations in the 1990s (through 1997) showed a surprising increase in both temperature and 
pressure (Elliot et ah, 1998). The next occultation was observed in 2008 but the data have yet 
to be reduced and analyzed, so the current state of Triton’s atmosphere is unknown (Sicardy 
et ah, 2008). Triton’s atmosphere is dominated by N 2 with traces of CO and CH4 (Tyler et 
ah, 1989; Lellouch et ah, 2010). Volatile transport, driven by migration of the sub-solar point, 
is responsible for the observed spatial distributions and temporal evolution of surface ices on 
Triton (Buratti et ah, 1994; Bauer et ah, 2010). Previously published papers describe in more 
detail the signihcant seasonal variations on Triton due to migration of the sub-solar point (e.g., 
Trafton, 1984; Hansen and Paige, 1992; Moore and Spencer, 1990). Fig. 1 shows the change in 
the sub-solar latitude over the period 1000-3000 C.E. A subset of that figure covering 1980-2020 
C.E. is presented in Fig. 2. After reaching its maximum southern extent in 2000 (-50°), the 
sub-solar pointj^has turned northward, reaching -42° at the time of this writing in mid-2015. 

* Triton is in a synchronous, retrograde orbit about Neptune. Its rotation is also retrograde with sunrise in 
the west. The sub-Neptune point is at 0° longitude. The south pole is the current summer pole. 
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Images obtained by Voyager 2 were taken during this extended period of southern illumination. 
A higher albedo in the southern hemisphere led some to argue for a south polar region covered 
in volatile ices (e.g., Stone and Miner, 1989; Moore and Spencer, 1990). With limited spectral 
information in the near-infrared from Voyager 2, the composition of various regions of Triton’s 
surface from the flyby epoch are unknown. 

As measured by Voyager 2, the mixing ratio of CH4 to N2 in Triton’s atmosphere is on the 
order of ~10“"^ (Broadfoot et ah, 1989; Herbert and Sandel, 1991). The presence of CH4, even 
at these levels, drives ongoing photochemistry (Lara et ah, 1997; Krasnopolsky and Cruikshank, 
1995). Indeed, photochemical haze was seen in Triton’s atmosphere by Voyager 2 (Herbert and 
Sandel, 1991; Rages and Pollack, 1992). The most common photochemical products of the 
interaction between CH4 and extreme-UV photons, cosmic rays, and charged particles from 
Neptune’s magnetosphere are acetylene (C2H2), ethylene (C2H4), and ethane (Krasnopolsky 
and Cruikshank, 1995; Moore and Hudson, 2003). Higher order hydrocarbons may also be 
created in much smaller quantities. Photochemistry occurs primarily in Triton’s atmosphere, 
as demonstrated by a calculation of the flux of Lyman-o; photons reaching the surface. The 
optical depth for Lyman-a photons in Triton’s atmosphere, r, is the product of the UV cross- 
section of CH4 at 120 nm (1.8x10“^'^ cm^; Chen and Wu, 2004) and the column density 
of CH4 (2.15x 10 ^®cm“^; Lellouch et al., 2010). This results in r=38.7 and means that the 
flux of Lyman-o; photons reaching Triton’s surface is so exceedingly small as to be negligible. 
Photochemical reactions involving CH4 occur exculsively in Triton’s atmosphere. 

Ethane, a photochemical product produced in large quantities (Krasnopolsky and Cruik¬ 
shank, 1995), was cited as a possible constituent ice on Triton’s surface by DeMeo et al. (2010). 
Theoretically, after forming in the atmosphere, non-volatile ethane would precipitate onto Tri¬ 
ton’s surface and accumulate over time. Once a 1 /rm spherical grain forms, it precipitates out 
of the atmosphere at the Stokes’ velocity (Mark Bullock, private communication): 

_ 2^/2 (p Patm) /.. \ 

where g is the gravitational acceleration on Triton (77.9 cm s“^), R is the radius of the ethane 
grains, p is the density of ethane (0.713 g cm“^; Barth and Toon, 2003), patm is the density 
of the atmosphere, and p is the dynamic viscosity of the atmosphere. The number density 
of the dominant component of Triton’s atmosphere, N 2 , is 10^® cm“^ at the surface (Herbert 
and Sandel, 1991), and we used this to calculate a mass density of 4.65x10“^ g cm“^. This 
value is negligible compared to the density of solid ethane, so it was ignored in the calculation. 
We calculate the dynamic viscosity using Sutherland’s formula (Eq. (2)). This gives a lower 
limit on the viscosity and therefore an upper limit to the Stokes’ velocity since the effects of 
turbulence and eddies are ignored. The dynamic viscosity (/i) of an N 2 atmosphere at T=50 K 
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Figure 1: Sub-solar latitude on Triton between 1000 and 3000 C.E. (data from JPL HORIZONS). The combi¬ 
nation of Neptune’s obliquity (30°), the inclination of Triton’s orbit (20°), and the rapid precession of Triton’s 
orbital node (637±40 years) contribute to Triton’s unique seasons (Trafton, 1984). This results in a beat pattern 
between the precession period and the 165-year orbit of Neptune. The shaded region covers the years 1980-2020 
and is presented in more detail in Fig. 2. 


(temperature of Triton’s troposphere; Tyler et al., 1989) is: 

To + C'/T^/^ 


( 2 ) 


where Tq is a reference temperature (300.55 K for N 2 ), /to is the dynamic viscosity at the 
reference temperature (0.0002 g cm“^ s“^ for N 2 ), and C is Sutherland’s constant for the gas in 
qnestion (111 K for N 2 ). We calcnlate /i=0.00003 g cm“^ s“^. Combining all this information 
yielded a descent velocity of ~0.004 cm s“^ for a 1 /tm ethane grain. Ethane more than 
likely does not form perfectly spherical grains, bnt instead takes the form of fractal aggregates; 
we assnme spheres for ease of calculation, so the descent velocity qnoted is an upper limit. 
Horizontal winds at altitndes below 8 km on Triton are 500 cm s ^ (Ingersoll, 1990). This 
means that within one equatorial circumnavigation of Triton, an ethane grain would descend, 
at most, on the order of 100 m from its altitnde of origin. Based on this calculation, we would 
expect ethane to precipitate nniformly onto Triton’s surface independent of where it formed in 
the atmosphere. With an accnmulation rate of 28 g cm“^ Gyr“^ (Krasnopolsky and Crnikshank, 
1995), we calculate that a total of 180 cm of ethane has precipitated uniformly onto the snrface 
of Triton over the age of the solar system. This assnmes that Triton’s atmosphere was in a 
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Figure 2: Sub-solar latitude on Triton between 1980 and 2020 C.E. (data from JPL HORIZONS). The maximum 
southern excursion of the sub-solar latitude was approximately -50° and occurred in 2000. The Voyager 2 
flyby took place when the sub-solar latitude was -45.43°. The lighter region denotes the time baseline of the 
IRTF/SpeX observations analyzed in this paper. The darker region denotes the range of sub-solar latitudes 
observed. The sub-solar latitude has monotonically increased over the period of our observations. 


steady state for the past 4.56 Gyr, so the calculated value is a rough upper limit. 

Ethane ice was previously detected on Pluto, an icy, outer solar system body similar in size 
to Triton (DeMeo et al., 2010; Holler et al., 2014). DeMeo et ah (2010) also made a tentative 
detection of ethane on Triton. In this work we study the longitudinal variability of the 2.405 /rm 
absorption feature and test whether or ethane is responsible for its presence. The 2.405 

/im feature is interesting because it could be due to ethane, ^^CO, or both (Cruikshank et ah, 
2006). DeMeo et al. (2010) hnd that the 2.405 pm feature in Triton’s spectrum is too strong 
to be due to ethane absorption alone; is therefore partially responsible for the observed 

depth. They also suggest the same for Pluto, presenting the possibility of a systematic error 
such as inaccurate optical constants. Work by Cruikshank et al. (2006) suggests that the 2.405 
pm band in the spectra of Pluto and Triton is due to ethane and not ^^CO. 

Sufficient solar illumination and the presence of strong vibrational transitions of N 2 , CH4, 
CO, CO 2 , and H 2 O make the near-infrared (0.7-2.5 pm) ideal for determining longitudinal 
distributions and for quantifying temporal evolution of these ices on the surface of Triton. 
We obtained reflectance spectra of Triton between 2002-2014 using the same telescope and 
instrument, allowing for easy comparison and combination of data in an effort to describe the 
effects of volatile transport on Triton’s surface in both space and time. We also investigated 
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the spatial and temporal variability of the absorption feature at 2.405 pm using this large data 
set. 

2. Observations &; Reduction 

Triton was observed from 2002-2014 using the SpeX infrared spectrograph at NASA’s 
3-meter Infrared Telescope Facility (IRTF) atop Manna Kea (Rayner et ah, 1998, 2003). 
During this period, 63 nights of usable reflectance spectra were obtained with observational 
circumstances found in Table 1 (information for 2002/07/17 to 2009/07/21 adapted from 
Grundy et ah, 2010). Text hies including the spectra for all good nights can be found at 
www2.lowell.edu/users/grundy/abstracts/2010.Triton.html and as supplementary material on¬ 
line. Spectra were obtained in the short cross-dispersed mode using a 0.3” xl5” slit. We aimed 
for a rotation angle close to the paralactic angle while avoiding angles within 15° of the imag¬ 
inary line connecting Neptune and Triton. Following this procedure signihcantly reduced hux 
contamination from Neptune. All spectra obtained prior to 2014 were resampled onto a consis¬ 
tent wavelength grid covering 0.67-2.55 pm at A/AA~2000 to match observations made with 
SpeX following its 2014 upgrade. 

Triton spectra were obtained using an ABBA slit dither pattern, with 2-minute integra¬ 
tions at each A and B position. Observations of Triton were interspersed with observations 
of a nearby solar analog star roughly once an hour. Solar analog spectra were also obtained 
following an ABBA dither pattern, with integration times dependent on the star. A solar 
analog requiring a large slew away from Neptune’s position would potentially result in a less- 
than-optimal airmass correction and cause instrumental hexure, so three different stars were 
used between 2002 and 2014. New solar analogs were found when Neptune’s motion across 
the sky produced too large a separation from the previous solar analog. From 2002-2007, HD 
202282 (spectral type G3 V; Honk and Smith-Moore, 1988) was used as the solar analog. From 
2007-2014, BS 8283 (binary with spectral types GO V and G1 IV; Neckel, 1986; Pourbaix et 
ah, 2004) was used as the solar analog. This binary star was checked against HD 202282 in 
2007 to ensure that it could serve as the solar analog. Starting in 2014, HD 215295 (spectral 
type G2/3 V; Honk and Swift, 1999) was used as the solar analog. 

An optimal extraction algorithm was used to reduce the Triton spectra (Horne, 1986). 
Wavelength calibration was performed using spectra of argon lamps contained within the SpeX 
instrument as well as sky emission lines in the Triton spectra themselves. The use of AB pairs 
eliminated sky emission lines. Airmass-corrected solar analog spectra were matched to Triton 
spectra with similar airmasses from the same night. Dividing the Triton spectra by the match¬ 
ing solar analog spectra removed solar absorption lines and corrected for telluric absorption. 
Some telluric absorption persisted because it is not possible for the Triton and solar analog 
spectra to be obtained simultaneously at the same airmass. All spectra were then normalized. 
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A grand average spectrum of all 63 individual nightly spectra was constructed by performing a 
weighted average within each wavelength bin (Fig. 3). Geometric albedos were not obtainable 
because the slit did not admit all light from the target object. 
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Figure 3: Combined Triton grand average spectrum of the 63 individual spectra obtained between 2002 and 
2014. Regions affected by significant telluric absorption were removed. The SpeX instrument currently covers 
0.67-2.55 /im at a spectral resolution of ~2000. All 63 spectra were resampled onto a consistent wavelength 
grid before being combined. The absorption bands considered in the analysis are identified and marked with 
dotted lines. 
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Table 1 

Observational Circumstances 


UT date 
mean-time 

Weather and 
i7-band seeing 

Sub-Earth 
Lon. (°E) 

Sub-Earth 
Lat. (°S) 

Phase 
angle (°) 

Total exp. 
time (min) 

2002/07/17 09:53 

Cirrus, 0.8” 

124.3 

49.8 

0.52 

84 

2002/07/18 09:36 

Thin cirrus, 0.8” 

184.8 

49.8 

0.49 

64 

2002/07/19 09:31 

Thin cirrus, 0.7” 

245.8 

49.8 

0.46 

68 

2002/07/20 09:27 

Thin cirrus, 0.8” 

306.9 

49.8 

0.42 

68 

2002/07/21 09:29 

Thin cirrus, 0.6” 

8.2 

49.8 

0.39 

68 

2002/07/22 09:21 

Clear, 0.7” 

69.1 

49.8 

0.36 

84 

2002/09/16 06:19 

Cirrus, 0.5” 

249.5 

49.9 

1.34 

84 

2002/10/03 06:32 

N/A, 0.6” 

211.0 

49.9 

1.67 

48 

2003/07/04 13:58 

Clear, 0.5” 

101.0 

49.6 

0.98 

64 

2003/08/10 10:09 

Some clouds, 0.5” 

196.4 

49.7 

0.19 

58 

2003/09/10 09:32 

Clear, 0.6” 

292.6 

49.8 

1.13 

40 

2003/10/16 05:51 

Patchy clouds, 0.5” 

327.5 

49.8 

1.80 

88 

2004/06/28 13:07 

Cirrus, 0.5” 

194.9 

49.3 

1.19 

152 

2004/08/12 12:19 

Clear, 0.7” 

67.7 

49.5 

0.21 

52 

2004/09/12 10:51 

Clear, 1.0” 

161.7 

49.6 

1.15 

48 

2004/10/21 06:13 

Partly cloudy, 0.8” 

18.0 

49.6 

1.84 

52 

2005/07/04 13:06 

Clear, 0.7” 

244.6 

48.9 

1.10 

146 

2005/09/19 08:02 

Clear, 0.5” 

265.2 

49.3 

1.26 

68 

2006/05/26 14:52 

Clear, 0.4” 

62.3 

48.4 

1.86 

40 

2006/06/26 13:39 

Cirrus, 0.4” 

157.9 

48.5 

1.37 

128 

2006/07/26 11:01 

Clear, 0.9” 

188.0 

48.6 

0.53 

72 

2006/08/30 09:14 

Clear, 0.7” 

166.0 

48.8 

0.62 

48 

2006/10/28 05:59 

Thin cirrus, 0.6” 

170.3 

49.0 

1.85 

68 

2007/06/21 13:15 

Clear, 0.8” 

252.6 

47.9 

1.53 

136 

2007/06/22 13:10 

Thin cirrus, 0.6” 

313.7 

47.9 

1.51 

160 

2007/06/23 13:18 

Clear, 0.6” 

15.2 

47.9 

1.49 

128 

2007/06/24 14:06 

Clear, 0.7” 

78.5 

47.9 

1.47 

88 

2007/06/25 13:13 

Clear, 0.5” 

137.5 

47.9 

1.45 

132 

2008/07/14 10:57 

Clear, 0.5” 

318.3 

47.4 

1.01 

72 

2008/09/01 07:02 

Clear, 0.6” 

68.1 

47.8 

0.55 

56 

2008/09/28 09:55 

Partly cloudy, 0.5” 

288.5 

48.0 

1.32 

60 

2009/06/22 14:35 

Clear, 0.7” 

102.2 

46.6 

1.59 

32 

2009/06/24 14:25 

Clear, 0.6” 

224.3 

46.6 

1.56 

80 

2009/07/21 12:07 

Clear, 0.5” 

71.7 

46.8 

0.88 

64 

2009/08/13 12:07 

Mostly clear, 0.6” 

39.8 

47.0 

0.15 

56 

2009/08/31 11:07 

N/A, 0.5” 

60.0 

47.2 

0.44 

42 

2010/08/07 11:08 

N/A, 0.7” 

71.9 

46.3 

0.43 

56 

2010/08/30 09:55 

Clear, 0.9” 

36.9 

46.5 

0.33 

80 

2011/06/11 15:01 

Clear, 0.5” 

233.0 

45.1 

1.83 

28 

2011/06/12 14:23 

Cirrus, 0.4” 

292.7 

45.1 

1.82 

86 

2011/06/20 14:40 

Clear, 1.0” 

63.4 

45.1 

1.72 

50 

2011/06/28 14:31 

Cirrus, 0.7” 

193.0 

45.2 

1.58 

60 

2011/08/06 12:04 

N/A, 0.6” 

54.8 

45.5 

0.55 

60 



Table 1 (cont.) 


Observational Circumstances 


UT date 
mean-time 

Weather and 
i7-band seeing 

Sub-Earth 
Lon. (°E) 

Sub-Earth 
Lat. (°S) 

Phase 
angle (°) 

Total exp. 
time (min) 

2011/08/23 09:08 

N/A, 0.5” 

8.1 

45.6 

0.02 

64 

2011/09/18 09:42 

N/A, 0.5” 

161.3 

45.9 

0.84 

56 

2011/11/13 06:07 

Mostly clear, 0.9” 

341.4 

46.2 

1.87 

144 

2012/06/25 14:19 

Clear, 0.7” 

111.6 

44.3 

1.66 

80 

2012/07/18 14:43 

Clear, 0.7” 

81.1 

44.5 

1.15 

66 

2012/09/02 07:01 

N/A, 0.7” 

357.9 

44.9 

0.29 

72 

2012/09/05 09:44 

Cirrus, 0.7” 

188.5 

45.0 

0.39 

284 

2012/09/08 08:41 

Cirrus, 0.7” 

9.5 

45.0 

0.49 

200 

2012/09/09 09:52 

Clear, 0.7” 

73.7 

45.0 

0.52 

224 

2012/09/10 09:51 

Clear, 0.6” 

134.9 

45.0 

0.55 

228 

2012/09/12 10:03 

Clear, 0.7” 

257.9 

45.1 

0.61 

232 

2012/09/13 08:54 

Clear, 0.7” 

316.1 

45.1 

0.64 

248 

2013/07/23 15:11 

N/A, 0.7” 

70.1 

43.6 

1.09 

36 

2013/07/25 14:55 

Clear, 0.9” 

191.9 

43.6 

1.03 

48 

2013/07/28 14:57 

N/A, 0.6” 

15.7 

43.6 

0.95 

56 

2014/08/04 13:49 

Mostly clear, 1.0” 

122.9 

42.8 

0.82 

135 

2014/08/05 13:19 

Clear, 0.7” 

182.9 

42.8 

0.79 

155 

2014/08/06 13:10 

Clear, 0.5” 

243.7 

42.8 

0.76 

175 

2014/08/10 12:57 

Clear, 0.7” 

128.1 

42.8 

0.63 

167 

2014/08/11 12:59 

Clear, 0.3” 

189.4 

42.9 

0.60 

171 


3. Analysis 

3.1. N 2 , CO, CH 4 , CO 2 , and the 2.405 /am band 

We calculated the integrated band areas (identical to an equivalent width calculation for a 
normalized spectrum) of the absorption bands presented in Table 2 in order to quantify N 2 , CO, 
CH 4 , CO 2 , and 2.405 /rm band variations with longitude and time on the surface of Triton. As 
a hrst step, we determined the wavelength range of three separate regions around each feature: 
pre-absorption continuum, band, and post-absorption continuum. These wavelength ranges are 
presented for each absorption band in Table 2 . A least squares linear fit was performed across 
the band on the points within the pre- and post-continuum regions (in this context, continuum 
does not refer to a region completely devoid of absorption). After the spectrum was normalized 
by the best-fit line, the integrated band area was calculated via numerical integration. The 
errors on the integrated band areas were computed in an identical manner to those in Grundy 
et ah ( 2010 ). 

The integrated band area of each night-averaged spectrum was plotted against its corre¬ 
sponding sub-observer longitude. A sinusoid of the form y = Asin(a; + B) + C, where y is 
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integrated band area (in /im) and x is longitude, was robustly fit to the longitudinal (zonal) 
distribution using a non-linear least squares algorithm. The data were initially £t to a sinusoid 
and the residuals and standard deviation of the residuals were calculated; data points with 
residuals greater than S-a were removed and another sinusoidal fit performed (Buie and Bus, 
1992). This robust-fitting process was repeated a total of three times. The period was fixed 
at 2n, representing one rotation of Triton. The reduced was calculated for the sinusoidal 
fit and a horizontal line placed at the mean value of C (to model zero variability over one 
Triton rotation). The number of degrees of freedom was 60 for the sinusoidal fit and 62 for 
the no-variability model. The longitude of peak absorption, peak-to-peak amplitude (twice the 
amplitude divided by the minimum), mean (C), and reduced values of the sinusoidal and 
no-variability models for each absorption band are presented in Table 3. 

To quantify temporal change, the integrated band area for each night-averaged spectrum 
was plotted against the corresponding Julian date. The observations from each night were 
made over a small range of sub-observer longitudes, so to disentangle the effects of longitudinal 
and temporal variability the value of the sinusoidal fit at the mean longitude for each night was 
subtracted off. The mean was added back on to maintain a similar mean value as the longitu¬ 
dinal fit. We plotted these values against phase angle for each absorption band to determine if 
the resulting distributions had a slope or followed a non-linear functional form. When plotted 
against phase angle, all distributions were essentially reasonable scatter about a mean, so phase 
angle effects were not considered any further in the analysis. A robust linear least squares fit 
was performed to determine if a particular ice species saw an increase or decrease in absorption 
between 2002 and 2014. The slope (m), slope detection level {m/am)) and y-intercept ( 6 ) for 
each absorption band are found in Table 4. 

We also investigated how combinations of CH 4 bands varied with longitude and time. In 
particular, we combined the weak CH4 bands (0.89, 1.15, and 1.33 /xm), strong CH4 bands 
(1.65, 1.73, and 2.20 /xm), and all CH4 bands. Combination involved addition of the integrated 
band areas for the appropriate CH4 bands on each night and propagation of errors. A com¬ 
bination of the bands comprising the CO 2 triplet (1.965, 2.01, and 2.07 /xm) was constructed 
in an identical manner. A combination of the 1.58 and 2.35 /rm CO bands was performed but 
the results did not add any new information. Results for the sinusoidal £t to longitude and the 
linear fit to time for each combination are found in Tables 3 and 4, respectively. 

3.2. H 2 O 

Analysis of the longitudinal distribution and temporal evolution of H 2 O ice on Triton pro¬ 
ceeded differently than for the other ices. As seen in Fig. 3, absorption due to H 2 O is very 
broad and overlaps with absorption bands of CH4, CO, and CO2. This prevented the calcula¬ 
tion of an integrated band area, so instead we calculated fractional band depth for H 2 O. The 
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fractional band depth was computed for each individual spectrum by dividing the mean value 
of the points in a relatively flat region of H 2 O absorption (1.5099-1.5694 pm) by the mean value 
of the points in a suitable continuum region (1.2195-1.2693 pm), and subtracting this quantity 
from one. While the fractional band depth is an even cruder estimation of abundance than the 
integrated band area, it still provides information on how H 2 O abundance differs with longitude 
and time. We aimed to determine differences and not absolute amounts of any species in this 
paper. The procedures described in the previous section were used to determine the longitudi¬ 
nal distribution and temporal evolution of H 2 O ice on the surface of Triton. Information on the 
sinusoidal £t to longitude and the linear £t to time are found in Tables 3 and 4, respectively. 

Table 2 


Continuum and Band Regions 


Absorption band 

Pre-absorption (pm) 

Band (pm) 

Post-absorption (pm) 

N 2 (2.15 pm) 

2.0925-2.1164 

2.1164-2.1739 

2.1739-2.1838 

CO (1.58 pm) 

1.5510-1.5701 

1.5701-1.5859 

1.5859-1.6047 

CO (2.35 pm) 

2.3336-2.3453 

2.3453-2.3592 

2.3592-2.3678 

CO 2 (1.965 pm) 

1.9564-1.9608 

1.9608-1.9707 

1.9707-1.9850 

CO 2 ( 2.01 pm) 

1.9995-2.0068 

2.0068-2.0178 

2.0178-2.0381 

CO 2 (2.07 pm) 

2.0539-2.0651 

2.0651-2.0736 

2.0736-2.0963 

2.405 pm 

2.3970-2.4014 

2.4014-2.4090 

2.4090-2.4134 

Weak CH4 bands: 

CH 4 (0.89 pm) 

0.8595-0.8776 

0.8776-0.9085 

0.9085-0.9298 

CH 4 (1.15 pm) 

1.0694-1.0995 

1.0995-1.2096 

1.2096-1.2296 

CH 4 (1.33 pm) 

1.2402-1.2920 

1.2920-1.3423 

1.3423-1.3435 

Strong CH4 bands: 

CH 4 (1.65 pm) 

1.6040-1.6194 

1.6194-1.6839 

1.6839-1.6985 

CH 4 (1.73 pm) 

1.6855-1.6970 

1.6970-1.7368 

1.7368-1.7439 

CH 4 (2.20 pm) 

2.1739-2.1838 

2.1838-2.2279 

2.2279-2.2360 
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Table 3 

Rotational Sinusoidal Fits 


Absorption 

band 

Peak 
long. (°) 

Peak-to-peak 
amp. (%) 

Mean 
( 10 “® /im) 

Reduced 

Sine Line 

N 2 (2.15 nm) 

35±1 

68±1 

3371.6±8.8 

2.8 

45 

CO (1.58 fim) 

60±3 

33±1 

342.6A1.5 

1.3 

2.1 

CO (2.35 fim) 

55±4 

90±6 

1016±13 

0.9 

3.7 

CO 2 (1.965 fim) 

97±54 

8±8 

732±20 

0.1 

0.1 

CO 2 ( 2.01 ytim) 

124±83 

4ih7 

846±20 

0.3 

0.3 

CO 2 (2.07 fxm) 

236A12 

15±3 

297.6A2.9 

0.2 

0.2 

CO 2 (All) 

128±61 

4±5 

1877A31 

0.3 

0.3 

H 2 OI] 

228A100 

1±1 

0.15938±0.00085 

1.5 

1.5 

2.405 fim 

152±44 

30±25 

95.3±7.1 

1.1 

1.1 

CH 4 (All Avg.) 

296±1 

22±1 

24144.8±3.2 

6.3 

48 

Weak CH 4 bands: 

CH 4 (Weak Avg.) 

285±1 

36±1 

8620.4±1.7 

15 

58 

CH 4 (0.89 /iui) 

304±1 

37±1 

566.3A1.7 

3.2 

5.0 

CH 4 (1.15 /iui) 

284±1 

37±1 

6140.2A4.3 

19 

54 

CH 4 (1.33 /ini) 

284±2 

32±1 

1901.9±5.6 

5.4 

14 

Strong CH 4 bands: 

CH 4 (Strong Avg.) 

307±2 

16±1 

15514±28 

3.2 

16.6 

CH 4 (1.65 /ini) 

313±2 

15±1 

8068±14 

2.4 

14 

CH 4 (1.73 /ini) 

315±4 

19±1 

3050±13 

1.4 

4.5 

CH 4 ( 2.20 /iui) 

289±1 

17±1 

4397.70A0.60 

1.9 

4.3 


tThe fractional band depth was calculated for the H 2 O absorption instead of the integrated band area. 
Therefore the mean value and uncertainty on the mean value are reported as unitless fractions. 
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Table 4 


Secular Linear Fits 


Absorption 

band 

m±am 
( 10 “® /im yr“^) 

m/Orn 

6 ± cTfe 

( 10 “® /im) 

Reduced 
m >0 m =0 

N 2 (2.15 fim) 

11.4±3.2 

3.6 

3373T12 

2.6 

2.8 

CO (1.58 /im) 

4.6±1.2 

4.0 

342.0±4.8 

1.0 

1.3 

CO (2.35 /im) 

-1.4±4.5 

0.3 

1015T18 

0.9 

0.9 

CO 2 (1.965 /im) 

2.6±4.6 

0.6 

732T19 

0.1 

0.1 

CO 2 ( 2.01 /im) 

6.1±5.0 

1.2 

848±21 

0.3 

0.3 

CO 2 (2.07 /im) 

1.0±2.4 

0.4 

297.5±9.7 

0.2 

0.2 

CO 2 (All) 

9.6±7.2 

1.3 

1879±30 

2.0 

2.1 

H20g 

O.OiO.O 

0.2 

0.15938±0.00084 

1.5 

1.5 

2.405 /im 

-2.4±1.8 

1.4 

94.4±7.3 

1.1 

1.1 

CH 4 (All) 

30.3±7.9 

3.8 

24154±32 

6.0 

6.1 

Weak CH4 bands: 

CH 4 (Weak) 

-32.1±4.3 

7.5 

8607T18 

14 

15 

CH 4 (0.89 /im) 

-3.4±1.3 

2.6 

566.1±5.5 

3.0 

3.3 

CH 4 (1.15 /im) 

-25.4±3.5 

7.3 

6131T14 

18 

19 

CH 4 (1.33 /im) 

-3.7±2.0 

1.9 

1899.5±8.0 

5.3 

5.3 

Strong CH4 bands: 

CH 4 (Strong) 

63.7±6.6 

9.7 

15529±27 

1.7 

3.1 

CH 4 (1.65 /im) 

23.1±3.4 

6.8 

8072T14 

1.6 

2.3 

CH 4 (1.73 /im) 

14.8±3.1 

4.7 

3053±23 

1.0 

1.4 

CH 4 ( 2.20 /im) 

26.4±4.7 

5.6 

4407T19 

1.3 

1.8 


4. Results Discussion 

The longitudinal (zonal) distributions of N2, CO, CH4 (all bands), CO2 (all bands), H2O, 
and the 2.405 /rm band are presented in Fig. 4-7 and 9-11. A discussion of temporal variability 
is found later in this section. Plots presenting absorption as a function of time are provided for 
CH4 (strong), N2, CO (2.35 /rm), and the CO2 triplet in Fig. 13-16. 

The two species with the largest peak-to-peak amplitudes are the volatiles N 2 (Fig. 4) and 
CO (1.58 /rm. Fig. 5; 2.35 fim, Fig. 6 ), with values of 68±1%, 33±1%, and 90±6%, respec¬ 
tively. The longitudes of peak absorption are 35±1°, 60±1°, and 55±4°, respectively. The 
phase difference between N 2 and the 1.58 /im CO band is 25±1°; between N 2 and the 2.35 fim 
CO band it is 20±4°; and between the two CO bands it is 5±5°. The longitudinal distributions 
of N 2 and CO are nearly in-phase, suggesting these species migrate together across the surface 


^The fractional band depth was calculated for the H 2 O absorption instead of the integrated band area. The 
slope and the uncertainty on the slope are instead reported in units of fractional band depth per year. The 
j/-intercept and the uncertainty on the j/-intercept are instead reported as unitless fractions. 
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of Triton. The large peak-to-peak amplitude as well as the relative and absolute longitudes of 
peak absorption of these species agree with previous work performed on a smaller SpeX data 
set in Grundy et ah (2010). Both ices have similar, non-negligible sublimation pressures (Fray 
and Schmitt, 2009) and are fully miscible in one another (Vetter et ah, 2007). However, note 
that the distributions are not completely in-phase, even when accounting for the uncertainty. 
The 1.58 /im CO band is found within a region of H2O absorption, but as described below, H2O 
ice shows negligible variability over one Triton rotation. The 2.35 /xm CO band is in between 
two strong CH4 bands at 2.32 and 2.38 /rm. However, the fact that the two CO bands are 
fully in-phase suggests that our choices of continua regions for each band eliminated any effect 
from the surrounding bands. Therefore, the peak longitudes of the CO bands are accurate and 
the shift with respect to N 2 is real. An explanation for this phenomenon is not immediately 
evident. 

Triton’s spectrum is dominated by CH4 absorption bands. We present here only the longi¬ 
tudinal variation of the sum of all CH4 bands (Fig. 7) since the individual CH4 bands all have 
roughly the same distribution. CH4 shows moderate variability over one Triton rotation with 
a peak-to-peak amplitude of 22±1%. The individual CH4 bands, and thus the groupings of 
bands, all have peak longitudes in the region between ~285-315° (Fig. 8), which is in agreement 
with the results from Grundy et ah (2010). These peak longitudes are found on the opposite 
side of the sub-Neptune hemisphere from those of N 2 and the two CO bands (phase differences 
of 99±1° from N 2 , 124±3° from the 1.58 /rm CO band, and 115±4° from the 2.35 yUm CO 
band). While only small amounts of CH4 can be diluted in N 2 , and vice versa (Prokhvatilov 
and Yantsevich, 1983), there is no reason to believe that the N 2 and CH4 distributions should 
be so far out of phase. In fact, it is possible for N 2 :CH 4 and CH 4 :N 2 mixtures to spatially 
coexist (Tegler et ah, 2012; Protopapa et ah, 2015). The fact that the two distributions are so 
far out of phase is unexpected. A similar pattern for the N 2 , CO, and CH4 distributions are 
seen on Pluto (Grundy et ah, 2013). This similarity may imply more than mere coincidence 
and requires further study. 

The non-volatile ices CO 2 and H 2 O show very little variation over one Triton rotation and 
have highly uncertain longitudes of peak absorption. We present the longitudinal variation of 
the sum of the CO 2 triplet (Fig. 9) rather than the individual CO 2 bands because the distribu¬ 
tions are similar and the sum has smaller error bars. The peak longitude for the combination of 
CO 2 bands is 128±61° and the peak-to-peak amplitude is 4±5%. We calculated the sinusoidal 
fit to the H 2 O fractional band depths (Fig. 10) and obtained a peak longitude of 228±100° 
and a peak-to-peak amplitude of 1±1%. The longitudinal distributions for both species can¬ 
not be distinguished between a sinusoid and a horizontal line (Table 3), meaning that their 
distributions across the surface of Triton are uniform. This is the same conclusion reached by 
Grundy et ah (2010) for the non-volatile ices. We note that the peak longitude (152±44°), low 
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amplitude, and reduced values for the 2.405 /rm band suggest a uniform surface distribution 
as well. 

We investigated what effects, if any, the two CH 4 absorption bands at 2.38 and 2.43 /rm had 
on the the longitudinal distribution of the 2.405 /rm band. The effect of CH 4 absorption must 
be negligible since the 2.405 /rm band distribution has such a low mean (0.0000953±0.0000071 
/rm) and a negligible amplitude (0.000012T0.000010 /rm). The peak-to-peak amplitude for the 

2.405 /xm band is potentially large (30±25%, rivaling the N 2 band), but this is misleading. The 
mean of the sinusoidal £t is very small and the amplitude is therefore large in comparison. 
However, the amplitude itself is low in an absolute sense. Contamination from more strongly 
variable CH 4 would be noticeable from its effects on both the mean and the amplitude of the 
sinusoidal fit. From this we conclude that our choice of continuum regions on either side of the 

2.405 /im band effectively removed contamination from the adjacent CH 4 bands. The species 
responsible for the 2.405 /rm band is distributed nearly nniformly across the surface of Triton. 

From the grand average spectrnm (Fig. 3) we detect the 2.405 /rm absorption band at the 
4.5-(T level. We compared the longitudinal distributions of CO (Fig. 5 & 6 ) and the 2.405 /rm 
band (Fig. 11) to determine the origin of the latter. The phase difference between the 2.405 
/xm band and the 1.58 /rm CO band is 92±44°. The phase difference between the 2.405 /rm 
band and the 2.35 /im CO band is 97±44°. Even if fractionation of CO were a factor on the 
snrface of Triton, we would not expect the longitudinal distributions of ^^CO and ^^CO to be 
this far ont of phase. Perhaps the quantity of ^^CO present on Triton is too low to detect 
in reflectance spectra. Assnming that carbon on Triton matches solar composition, 1% of all 
CO on Triton should be in the form of ^^CO (Scott et ah, 2006). Even if ^^CO absorption at 

2.405 /im were 1% as strong as ^^CO absorption at 2.35 /rm (Fig. 3), the feature would be 
indistinguishable from noise. The longitndinal distribntion of the 2.405 /im band also shares 
featnres with the distributions of non-volatile CO 2 and H 2 O (Table 3): a poorly constrained 
peak longitnde and low amplitude. This provides additional circumstantial evidence that the 
species responsible for the 2.405 /im band is also relatively non-volatile. Both ^^CO and ^^CO 
are volatile whereas ethane is non-volatile. This snpports the interpretation that ethane is 
responsible for the absorption feature at 2.405 /im. 

A summary of the longitudinal distributions across Triton’s snrface is presented in Fig. 12. 
From this figure we note the large nncertainties on the peak longitndes for the non-volatile ices 
(CO 2 , H 2 O, ethane) and the smaller uncertainties on the peak longitudes for the volatile ices 
(N 2 , CO, CH4). The volatile ices also show larger variability across the surface of Triton, with 
the longitnde of peak absorption likely corresponding to a large patch of a particular ice species. 
CO 2 and H 2 O are uniformly distribnted across the snrface and act as the snbstrate upon which 
the volatile ices non-uniformly deposit (with ethane as a contaminant). One would expect the 
deposition of volatile ices to indnce variability in the non-volatile ices, especially in regions 
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of thick volatile coverage, yet non-volatile variability is negligible. The viewing geometry of 
Triton provides an explanation: Any latitude less than 90° minus the sub-solar latitude will be 
visible at all times. The simplest explanation for the observed variability of non-volatile ices is 
that they dominate the always visible south polar region. Therefore, the south polar region is 
relatively devoid of volatile ices. The volatile ices should be found farther north where parts of 
the surface rotate into and out of view. 

We support this claim with our analysis of the temporal evolution of the volatile and non¬ 
volatile ices. Statistically significant slopes (>3-(j) were found for the strong CH 4 band group¬ 
ing, N2, CO (1.58 /im), and CH4 (all, weak band grouping, 1.15 /rm, 1.65 /rm, 1.73 /rm, 2.20 
fim). However, for all but the strong CH4 band grouping (Fig. 13), the reduced values 
do not favor a sloped line over a horizontal line (Table 4). We find that a positively sloped 
line is statistically preferred for the strong CH4 band grouping, meaning the quantity of CH4 
present on the visible portion of Triton increased over time. Figures 14-16 present the tem¬ 
poral distributions for N2, CO (2.35 /rm), and the CO2 triplet. N2 and CO are volatile ices 
that might be expected to show a similar change in time to volatile CH4; the two CO bands 
show similar changes in time, so only one is presented. From our analysis, we conclude that 
the signal-to-noise ratios of the nightly spectra are not high enough to distinguish between a 
flat and non-zero slope for any ices other than CH4. Grundy et ah (2010) found that the N2 
ice integrated band area showed a small decrease between 2000 and 2009 by considering two 
subsets of their data: 2000-2004 and 2005-2009 (see Fig. 11 of that paper). A slight downward 
shift in the longitudinal distribution was detected over this time period that they interpret as 
a change in the texture of the N 2 ice, not a change in areal coverage. A negligible change in 
N 2 absorption, as reported by Grundy et ah (2010), is consistent with the results of this work. 
Analysis of temporal variation for other species was not provided in their work. 

We considered the possibility that the observed increase in CH4 absorption was due to 
volatile transport. Evidence for volatile transport on Triton dates to the 1950s, with conclusive 
evidence for surface color change between 1979 and the Voyager 2 flyby in 1989 (Smith et ah, 
1989; Buratti et ah, 1994; Brown et ah, 1995). The disk-averaged visible spectrum (0.3-0 .8 fj,m) 
of Triton became significantly less red between 1979 and 1989 as the sub-solar latitude moved 
10° south. This color change was interpreted as deposition of volatile ices onto older, redder 
regions of the surface. Bauer et ah (2010) present evidence for significant albedo changes 
across the surface of Triton between the Voyager 2 flyby in 1989 and 2005. The albedo of the 
equatorial regions and the sub-Neptune hemisphere increased and the albedo decreased in the 
anti-Neptune hemisphere over this period. We note that volatile ice absorptions peak on the 
sub-Neptune hemisphere (Fig. 12, this work). Fig. 2 from Bauer et ah (2010) shows the 2005 
visible albedo map compared to the Voyager 2 albedo map. Shortly after the Voyager 2 flyby in 
1989, many argued that volatile ices were found in the southern hemisphere due to the higher 
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albedo (e.g., Stone and Miner, 1989; Moore and Spencer, 1990). This remains unconfirmed due 
to limited spectral information from Voyager 2. 

The increase in CH 4 absorption may instead be due to a change in viewing geometry. Our 
hndings suggest that the southern latitudes are primarily denuded of volatile ices, revealing the 
underlying non-volatile substrate. We propose that volatile ices migrated during the southern 
summer to the northern hemisphere where they deposited in the region near the summer/winter 
terminator. As the sub-solar (and thus, the sub-observer) latitude moves northwards, we re¬ 
ceive more direct reflection from these northern latitudes and less from the region surrounding 
the south pole. We dehne the south polar region to be the entire southern hemisphere due to 
its visual characteristics (Moore and Spencer, 1990). The northern temperate region is defined 
between 0° and the most northerly observable latitude. Imagine the sub-observer point at Tri¬ 
ton’s equator; the north pole and south pole are 90° away from the equator and just visible on 
the edge of the disk. For other values of the sub-observer latitude, the visible region will be 
bounded by the nearest pole at one extreme and the latitude 90° away on the other. 

At the start of the observations described in this work (July 2002 ), the sub-observer latitude 
was -49.86°; on the night of the last observation in August 2014, it was -43.15°. In both cases, 
the most southerly observable latitude was -90°. In 2002, the most northerly observable latitude 
was 40.14°; in 2014 it was 46.85°. We developed code that calculated the projected solid angle 
of each small area of Triton’s surface to determine the total projected area of a chosen range 
of latitudes. Then we divided this value by the total projected area (equal to tt) to obtain the 
percentage of Triton’s visible disk due to the chosen region. Due to longitudinal symmetry, this 
percentage applies at all rotational phases. See the Appendix for a more detailed mathematical 
description. The south polar region was responsible for 88.2% of the projected area in 2002 and 
84.2% in 2014, a fractional decrease of 0.05. The northern temperate region was responsible 
for 11.8% of the projected area in 2002 and 15.8% in 2014, a fractional increase of 0.34. 

A signihcant increase in albedo near the longitude of peak absorption for CH4 is seen in 
the 2005 map from Bauer et al. (2010). No such brightening is noticed at the longitudes of 
peak N 2 or CO absorption and the latitudes of increased albedo were visible in both 1989 and 
2005 (this region was visible for the entire duration of our observations). These facts cannot be 
explained by a change in viewing geometry and a volatile transport explanation is not imme¬ 
diately evident. A variety of future observations and modeling of volatile transport on Triton 
will be necessary to explain the observed temporal changes. 
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Figure 4: Integrated band area of the 2.15 /rm N 2 band plotted against sub-observer longitude. The dashed 
line is a robust sinusoidal fit to the points. N 2 ice shows significant variation over one Triton rotation. 



Figure 5: Integrated band area of the 1.58 /im CO band plotted against sub-observer longitude. The dashed 
line is a robust sinusoidal fit to the points. The phase difference between the 1.58 /rm and 2.35 fim CO bands 
is 5±5°. 
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Figure 6 : Integrated band area of the 2.35 /rm CO band plotted against sub-observer longitude. The dashed 
line is a robust sinusoidal fit to the points. CO ice shows significant variation over one Triton rotation with a 
peak longitude similar to that of N 2 ice. This suggests that N 2 and CO migrate together as Triton’s seasons 
change. The phase difference between the two CO bands is 5±5°. 



Figure 7: Total integrated band area of all CH 4 bands plotted against sub-observer longitude. The dashed line 
is a robust sinusoidal fit to the points. CH 4 ice shows moderate variation over one Triton rotation. 
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Figure 8 : Range of longitude of peak absorption (longitude ± uncertainty) for all CH 4 bands and combination 
of CH 4 bands. Bar sizes do not represent the extent of the ices but instead show the possible location of the 
peak longitude for that species. The 0.89, 1.15, and 1.33 /j,m bands comprise the weak CH 4 grouping. The 1.65, 
1.73, and 2.20 /rm bands comprise the strong CH 4 grouping. Triton is viewed looking down on the south pole, 
with rotation and orbital motion counterclockwise. 



East Longitude (°) 


Figure 9: Total integrated band area of all CO 2 bands plotted against sub-observer longitude. The dashed line 
is a robust sinusoidal fit to the points. CO 2 ice shows negligible variation over one Triton rotation. Additionally, 
the longitude of peak absorption for CO 2 is uncertain (Table 3). 


20 




































































Figure 10: Fractional band depth of H 2 O ice absorption plotted against sub-observer longitude. The dashed 
line is a robust sinusoidal fit to the points. H 2 O ice shows negligible variation over one Triton rotation with a 
large uncertainty in the longitude of peak absorption (Table 3). 



East Longitude (°) 

Figure 11: Integrated band area of the 2.405 /rm band plotted against sub-observer longitude. The dashed line 
is a robust sinusoidal fit to the points. The 2.405 p,m band shows negligible variation over one Triton rotation 
with large uncertainties in both the longitude of peak absorption and peak-to-peak amplitude (Table 3). The 
2.405 fim band is consistent with ethane absorption. Negative values of integrated band area correspond to 
individual nightly spectra with lower signal-to-noise ratios and have no physical meaning. 
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Figure 12: Range of longitude of peak absorption (longitude ± uncertainty) for all ice species found on Triton’s 
surface. The peaks from the longitudinal distributions of the combined CH 4 bands and the combined CO 2 
bands are presented for those species (denoted by “All”). Bar sizes do not represent the extent of the ices but 
instead show the possible location of the peak longitude for that species. Legends are in the quadrant of their 
respective distributions. Triton is viewed looking down on the south pole, with rotation and orbital motion 
counterclockwise. 
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Figure 13: Integrated band area of the combination of strong CH 4 bands plotted against Julian date. The 
dashed line is the robust linear fit to the points. We found an increase in absorption for the combination of 
strong CH 4 bands over the period of our observations at the 9.7-ct level. Reduced values (Table 4) favor a 
sloped line over a horizontal line. The integrated band area difference was calculated by subtracting the value 
of the sinusoidal fit at the proper phase from the integrated band area; the mean was added to maintain the 
same offset as the longitudinal distribution. 
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Figure 14: Integrated band area of the N 2 band plotted against Julian date. The dashed line is the robust 
linear fit to the points. We found a positive slope at the 3.6-tT level but the reduced values (Table 4) do not 
favor a sloped line over a horizontal line. The integrated band area difference was calculated by subtracting the 
value of the sinusoidal fit at the proper phase from the integrated band area; the mean was added to maintain 
the same offset as the longitudinal distribution. 
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Julian Date + 2452000 

Figure 15: Integrated band area of the 2.35 CO band plotted against Julian date. The dashed line is the 
robust linear fit to the points. We found a negative slope at the 0.3-cr level and the reduced values (Table 
4) do not favor a sloped line over a horizontal line. The integrated band area difference was calculated by 
subtracting the value of the sinusoidal fit at the proper phase from the integrated band area; the mean was 
added to maintain the same offset as the longitudinal distribution. 
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Figure 16: Integrated band area of the combination of CO 2 bands plotted against Julian date. The dashed line 
is the robust linear fit to the points. We found a positive slope at the 1.3 -ct level and the reduced values 
(Table 4) do not favor a sloped line over a horizontal line. The integrated band area difference was calculated 
by subtracting the value of the sinusoidal fit at the proper phase from the integrated band area; the mean was 
added to maintain the same offset as the longitudinal distribution. 
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5. Summary 

Between 2002 and 2014, 63 near-infrared spectra were obtained with the SpeX instrument 
at NASA’s IRTF. Analysis of the longitudinal (zonal) and temporal distributions of various 
ices indicate compositional changes in the 25 years since the Voyager 2 flyby. The hndings of 
this work are summarized below: 

• The volatile ices, N 2 , CO, and CH4, show signihcant variability over one Triton rotation. 
N 2 and CO have similar longitudes of peak absorption, with CH4 offset by about 90°. 

• The non-volatile ices, CO 2 and H 2 O, show little variability over one Triton rotation. 
In fact, the longitudinal distributions for these ices are consistent with horizontal lines, 
suggesting a uniform distribution across the visible surface of Triton. 

• The 2.405 /xm band was detected at the 4.5-a level in the grand average spectrum. Both 
non-volatile ethane (C 2 H 6 ) and volatile ^^CO absorb at 2.405 /rm. The longitudinal dis¬ 
tribution of the 2.405 /rm band is 92±44° and 97±44° out of phase with the distributions 
for the 1.58 /rm and 2.35 fim CO bands, respectively. If there is negligible fractionation 
of CO, this implies that the band cannot be due to ^^CO. Additionally, the longitudinal 
distribution of the band cannot be distinguished between a sinusoid and a horizontal line, 
suggesting a uniform distribution across the surface, similar to the other non-volatile ices 
CO 2 and H 2 O. This evidence is consistent with ethane absorption at 2.405 /rm. 

• Over the period of the observations, absorption from the combination of strong CH4 bands 
(1.65, 1.73, and 2.20 /rm) increased. While slopes were detected above the 3.0-(T level for 
N 2 , CO (1.58 /xm), and CH4 (all, weak, 1.15 /rm, 1.65 fj,m, 1.73 fj,m, 2.20 /rm), the reduced 

values for these bands do not favor a sloped line over a horizontal line. The increase 
in CH 4 absorption between 2002 and 2014 may be due to volatile transport or a change 
in viewing geometry. Slopes for all other absorption bands were detected below the 3.0-(j 
level and are therefore consistent with the null hypothesis of no temporal change. 

• The Voyager 2 flyby provided images of Triton’s southern latitudes, but no spectra. The 
high albedo of the southern hemisphere in 1989 suggested coverage by volatile ices. The 
combination of volatile and non-volatile ice distributions implies a southern hemisphere 
currently denuded of volatile ices. Revealing the full extent of temporal changes since the 
Voyager 2 flyby will be the subject of future work. 

IRTF/SpeX observations of Triton are ongoing, with hve nights of spectra obtained in 2015 
and a proposal submitted to continue the program in 2016. New observations are required to 
fully understand the processes that are altering Triton’s surface. Spatially-resolved spectra of 
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Triton would allow us to determine the latitudes of volatile and non-volatile ices. An albedo 
map of Triton, for comparison to previous maps made in 1989 and 2005, would be useful for 
quantifying the movement of ices across the surface. Information from occultations would 
indicate whether the surface pressure increased or decreased since the last observed occultation 
in 1997. Together, this information would allow us to understand Triton’s surface evolution 
since the Voyager 2 flyby in 1989 and to predict future changes. 
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Appendix 

We derive here an equation for the calculation of the projected solid angle of a portion of a 
sphere. An inhnitesimal solid angle (dG) on the surface of a sphere can be calculated from the 
product of the inhnitesimal latitude (d6*) and longitude (d0), and the cosine of the latitude: 

dG = d6'd0cos6*. (Al) 

The cos6 term is necessary since the solid angle of the poles themselves is zero. Integrating Eq. 
(Al) over the entire surface of the sphere (0 from 0 to 27r and 6 from — | to |) results in a 
solid angle of 47r, as expected. 

For an observer, the emission angle (y) from each dG affects the calculation of the surface 
area visible at any given time. The projected solid angle {flproj) is the integral of the product 
of the solid angle and the cosine of the emission angle: 

r02 j‘S2 

^proj = / / cosxcos6*d6'd0. (A2) 

J (Pi J 9i 

The total projected solid angle has a value of tt and can be visualized as an arbitrary hemisphere 
projected onto a circular disk of the same radius. 

Now consider two vectors, r and fj), that represent the direction from the center of the 
sphere to a surface element dG and the direction from the center of the sphere to the observer, 
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respectively. The x, y, and 2 : components of r and fo are measured from the center of the 
sphere as seen by the observer. In vector notation, we have: 


r = (cosOcoscj), cosOsmcj), sin 6 ') 


(A3) 


and 

vq = (cos6*ocos0o, cos6*osin(;/)o, sin6*o). 

The emission angle is dehned by the dot product of Eq. (A3) and (A4): 


(A4) 


r ■ fo = |r||fo|cosx cosy = 


r ■ ro 


(AS) 


Since r and Fq are unit vectors, their moduli are 1. Therefore, cosy is equal to the dot product 
of the two unit vectors: 


cosy = r ■ f)) = cos6*cos0cos6*ocos0o + cos 6 *sin(;/)cos 6 *osin(;/)o + sin 6 *sin 6 'o (AS) 

The sub-observer longitude is chosen to be zero by rotational symmetry, simplifying the ex¬ 
pression for cosy: 

cosy = cos6*cos0cos6'o -|- sin 6 *sin 6 *o (A 6 ) 

Plugging Eq. (A 6 ) into Eq. (A2) yields the projected area in terms of latitude, longitude, 
sub-observer latitude, and sub-observer longitude: 

r4>2 j'd2 

^proj = / / cos 6 '[cos 6 *cos())cos 6 'o -l- sin0sin6*o]d6'd(;/). (A7) 

J (pi J 9i 

In this work we implemented a numerical calculation of Eq. (A7) in IDL. We started by 
considered all longitudes (0i=O and <p2=27r) for a chosen range of latitudes, then ignored any 
inhnitesimal surface areas with a corresponding value of cosy that was negative. By dividing 
the resulting projected area by tt, the value of the total projected area, we obtained a percentage 
of the projected area due to the region in question. 
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